A renormalized excitonic method in terms of block excitations. Application to spin 

lattices 
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Dividing the lattice into blocks with singlet ground state and knowing the exact low energy 
spectrum of the blocks and of dimers (or trimers) of blocks, it is possible to approach the lowest 
part of the lattice spectrum through an excitonic type effective model. The potentialities of the 
method are illustrated on the 1-D frustrated chain and the 1/5-depleted square and the plaquette 
2-D lattices. The method correctly locates the phase transitions between gapped and non-gapped 
phases. 



I. INTRODUCTION 

The idea that one may see a periodic lattice as built 
from interacting blocks of sites, rather than as interacting 
sites, is computationally and intellectually attractive and 
has received several exploitations. This process consists 
in a scale change. Wilson's bright proposal of real space 
renormalization group (RSRG)i is certainly the most el- 
egant illustration of this idea, since it can be infinitely it- 
erated and asymptotically converges on physically mean- 
ingful accumulation points. The method in its original 
version consists in a simple and severe truncation of the 
Hilbert space by considering only the lowest states of each 
block and the products of the selected block eigenfunc- 
tions to approach the eigenfunctions of blocks of blocks. 
The method happened to be numerically inefficient. Re- 
cent works have shown that its performances can be dra- 
matically improved when one defines effective interac- 
tions between the blocks. These effective interactions are 
calculated from the knowledge of the exact spectrum of 
the dimers (or trimers) of blocks using Bloch's theory 2 
of effective Hamiltonian. The so-called contractor renor- 
malization (CORE) technique2*^&*£ considers blocks, re- 
tains a few eigenstates \Ia) of each block A and defines 
inter-block effective interactions (I aJ B\H e f ' \K aL b) , us- 
ing the spectrum of the AB problem and the Bloch's 
theory 2 of effective Hamiltonians. In most applications 
the size of the blocks remains small, several states per 
blocks are kept, and three and / or four blocks effec- 
tive interactions are introduced from the knowledge of 
the spectrum of trimers and / or tetramers of blocks. A 
specific variant of the method has been proposed by two 
of the authors under the name of RSRG-El£&& (RSRG 
with effective interactions) which treats spin lattices by 
considering blocks with an odd number of sites and a 
doublet ground state which is the only one to be explic- 
itily retained. The blocks are then quasi-spin. From the 
spectrum of dimers or trimers of blocks one may define 
an inter-block Heisenberg Hamiltonian. A proper design 
of the blocks frequently results in an isomorphism be- 
tween the original lattice and the lattice of blocks. Hence 
the process may be iterated, exhibiting critical ratios 
of the elementary interactions and accumulation points. 
The methodological studies have examined the following 



dilemma in the research of accuracy 

- consider larger blocks and only dimers of blocks 
(i.e., two-body effective interactions only) or 

- consider smaller blocks and trimers or tetramers 
(i.e., three and four blocks interactions) 

showing that in most cases the former solution is suf- 
ficient. In a rather similar spirit one may mention the 
block correlated coupled cluster method^ which employs 
the coupled cluster formalis m 1 ^ 12 ! 13 starting from the 
product of the ground state eigenfunctions for each block. 
A similar philosophy (with small blocks) is present in the 
applications of self-consistent perturbation formalism to 
periodic lattices^ These last two methods do not pro- 
vide information on the gaps, while CORE and RSRG-EI 
bring good estimates of them. 

The present work is closely related but is focused on a 
direct research of the gap. One starts now from blocks 
constituted of an even number of sites and having a non- 
degenerate singlet ground state. The blocks may be iden- 
tical or not, but they must lead to a periodic picture 
of the lattice in terms of blocks (hence with larger unit 
cells). The ground state will be built from products of 
block ground states. Considering the exact energies of 
dimers or trimers of blocks, one will define effective inter- 
actions between blocks in their ground states, producing 
an additive energy systematics. 

For the study of excited states, one will also con- 
sider the lowest states of the blocks and of the various 
dimers and eventually trimers of blocks. The knowledge 
of the lowest states of the dimers and trimers (energy 
and eigenvectors) enables one to define the effective in- 
teractions between an excited block and neighbor ground 
state blocks and effective excitation transfer integrals 
from one block to its neighbors. These quantities will 
be defined through the Bloch's theory of effective Hamil- 
tonians. Then the lattice lowest excitations are treated 
through an excitonic model that makes use of these ef- 
fective quantities. The theory is developed in section 
II. Section III first shows the improvement brought by 
the use of effective interactions rather than of bare inter- 
actions on elementary mono-electronic problems where 
the original RSRG version failed. The efficiency of the 
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method will then be illustrated on three spin lattices, 
namely the 1-D frustrated chain, the 1/5-depleted 2-D 
square lattice and the plaquette lattice. The three prob- 
lems exhibit phase transitions (of second and first order) 
which are satisfactorily treated with the here-proposed 
renormalized excitonic method (REM). 

II. METHOD 
A. Principle 

Let us consider a lattice constituted of blocks A,B... 
having a non-degenerate singlet ground state. %jP A is the 
ground state for the block A, of Hamiltonian Ha, 

H A iP° A = EM. (1) 

The zero-order description of the ground state of the lat- 
tice will be the product of the ground states for each 
block 

*o=n^- ( 2 ) 

A 

The zero-order energy would be additive 

£°=5>1, (3) 

A 

but since the Hamiltonian involves interaction operators 
between blocks 

ff = E^ + EE^ ( 4 ) 

A A <B 

the mean energy implies inter-block interactions 

A A <B 

= E £ a+EE»v (5) 

A A <B 

Here v AB is a zero-order interaction energy between 
blocks A and B. 

Solving exactly the AB problem 

{H A + H%+V AB )\^ AB )=E AB \* AB ), (6) 

enables one to define an improved interaction energy vab 

v AB =E AB -E A -E° B , (7) 

which takes into account at all orders the perturbative ef- 
fect of excitations on A and B, i.e., of the vectors \4> A 4> B )- 
Going to trimers it would be possible to define a quantity 

v AB c = E° ABC - E° A - E% - E% - v AB - v BC - v AC . (8) 



The ground state energy will be additive 

A A <B A <B <C 

At this stage we would like to stress on the fact that the 
treatment of the ground state is from a single function 
*o given by (Eq. 0. 

For the description of the excited states we shall con- 
sider a very limited model space. If one considers the 
lowest excited state on A, ^p* A , of the desired spin multi- 
plicity, 

HaWa) = E* A \p A ), (io) 
Wa) = T+|Vft>, 

the on-block excitation lowest energy is = E* A — E\. 
We intend to represent the lowest excitations on an en- 
semble of blocks from localy singly excited configurations 
of the type = T+* = \i/j a ■ ■ ■ tp° H ^i^°j ■ ■ ■ \ where the 
excitation is localized on block /. This will lead to an 
excitonic treatment of the excitation 




In order to establish the corresponding model Hamilto- 
nian one needs to calculate 

- the effective interaction between an excited block 
and neighbor blocks in their ground state 

- the excitation hopping integrals which propagate 
an excitation from one block to other blocks. 

These informations will be obtained from the spectral 
properties of dimers and trimers of blocks, using the ef- 
fective Hamiltonian theory. Let us consider first the ex- 
traction of information from interacting pairs of blocks. 

B. Extraction of information from dimers of blocks 

For a dimer AB one may chose as model space for the 
lowest energy excited states the two vectors tp A ip B and 
ip A ip B - The corresponding projector is 

p*ab = I^X^bI + I^IX^bI (12) 
= |r+* )(T+*°| + |r+*°)(T+^°|. 

If one identifies the two eigenvectors ^* AB and ^* AB of 
eigenenergies E* AB and E* AB 

Hab\^ ab ) = E*ab\**ab), (13) 
Habeas) = (14) 

which have the largest projections onto the model space, 
it is possible to define an effective Hamiltonian built on 
the model space and according to Bloch's definition 

H eff \pi B n B ) = E* AB \px B n B ), (is) 
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H eff \P%B*AB) = Eab\Pab*ab) 



(16) 



In order to have an hermitian effective Hamiltonian its 
eigenvectors must be orthogonal. We shall assume that 
is orthogonal or Schmidt-orthogonalized to 
One may write, after normalization, 



P*ab\^*ab 
P%b\**ab 



\Pab**ab) 



a\p A 1>%) 



\P1b*ab) = -Wa^b) + a\iP A r B ). 



(17) 



(18) 



The spectral definition of H e ^ leads to the following 
equations 



a'E AB -+ 
E* A + Ef 



b z E 



ab 

V( A * )B , (19) 



with equal or different respective interactions. The 
ground state being represented by in^l nas an ener gy 

K 



(28) 



K 



K <L 



The set of excited states described by our method will be 
built in the space spanned by all products 



if} = T/*° 



n n ^ 



. K =0,7-1 



,L=I+l, 



Their effective energy is 

£* = {^}\H eff \^}) = E°k+ e *i 



(29) 



E E ^ + E u ^*)- ( 3 °) 

K&I) <L(^I) K±I 



b E 
E° A 



ab 



a 2 E 



AD 



Et 



VA(B'), (20) 



The local excitation energy is 

e\ -s = m - ^? + E ( v «(i*) - v ki)- (si) 



(r A r B \H eff \r A rB) = (e 



AB 



E* AB )ab 



h 



AB- 



(21) 



The terms (resp. u^fs*)) represent the effective 

interactions between A* and B (resp. between A and 
B*) and h A B is the effective interaction responsible for 
the transfer of excitation from A to B. If A and B are 
identical blocks and if the AB dimer presents an element 
of symmetry transforming A in to B and vice versa, \a\ = 
\b\ = l/v2, one eigenvector is an in-phase combination 
of ^ a iPb an d V'aV's) °f energy E* AB , the other one being 
the out-of-phase combination, of energy E* AB . 



^AB^AB 



H AB \f> 



ABl 



1 



AB I ABl ' 



(22) 
(23) 



The vectors if!} and iffj interact through the matrix el- 
ement (T+ff°\H e ff\Tfff°) = (f}\H e ff\f*j) = hu. The 
effective Hamiltonian matrix has a near-diagonal struc- 
ture, similar to that of a tight-binding mono-electronic 
Hamiltonian. It generates bands which only represent 
the states of the lattice having large projections onto the 
vectors iff}, i.e., on the intra-blocks lowest energy excita- 
tions. The descriptions of the lowest energy states of the 
lattice should be relevant. If the blocks are identical and 
engaged in the same interactions of negative sign with 
their first neighbors, the excitation energy to the lowest 
k = state should be 



V K l) + 



h 



IK- 



(32) 



K=£I 



If the KI couples present an element of symmetry trans- 
forming K into /, using Eq. ESand one obtains, 



Pab^ab 



HA B \n u B ) = 



']{A»)B = VA(B') 



h A i 



V2 



E* A u B \n u B ) 



; (E A %+ E */ B 



)-E* 



E B , 



-(E A %~ET B ). 



(24) 
(25) 

(26) 

(27) 



It is then possible to consider the infinite lattice in which 
each block is surrounded by nearest-neighbor blocks B 



K 



(A 



KI 



(33) 



Ekj are excitation 



where A} = Ef - E°, A* KI = E*gj 
energies on the blocks and dimers of blocks respectively. 
One notices that the other root E^j of the dimer disap- 
pears in this expression. 

One sees that the derivation leads to a renormalized 
excitonic method, where the excitation transfer inte- 
grals hu are renormalized, therefore including to all or- 
ders some indirect processes going through higher-energy 
(multiple) excitations on neighbor blocks or inter-block 
excitations as will be shown hereafter. Of course the re- 
sults will be dependent on the shape and size n of the 
blocks. 
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C. Extraction of information from trimers of blocks 

It is possible to use the eigenstates of trimers of blocks 
to extract three blocks interactions. For a given shape 
of the elementary blocks one must of course consider the 
various types of trimers of blocks. While for the ground 
state the three-block correction is given by 

vabc = E abc -E a -E b -E%-v A b-v ac -vbc, (34) 

for the excited states the model space involves three vec- 
tors. The projector on the model space is 

+IVMV>c>(VM^I 

= |T+*°)(T+* | + |T+*°>(T+*°| 

+ |T+v&°)(T+vI< |, (35) 

one must identify the three eigenstates of the ABC prob- 
lem having the largest projections on the model space. 
Notice that these three states are not necessarily the 
three lowest ones. Let us call *& A bc ^*abc an d ^*abc 
these three states and E ABC , E* ABC and E* ABC the cor- 
responding eigenenergies. For hermiticity the projected 
eigenvectors are orthonormalized leading to three vec- 
tors $* ABC = PXbcVabo ^Ibc = P*abc^*abc and 



(T+y \H eff \T+V ) 



(T+V \H eff \T+y°) 



The excitation energy for the vector k = (which is not 
necessarily the lowest one), when all blocks are equiva- 
lent, is 

A*°° = E*j - Ej + 5>(/-)., - v u ) + 
j 

^2{V(I*)JK - vijk) + ^2(hu + h u(K))- (41) 

JK J K 

The method is generalizable to four (and more) blocks. 
One should however remark that when one increases the 
number of blocks the identification of the eigenstates hav- 
ing the largest projections onto the model space may be- 
come ambiguous. When changing the ratio of the inter- 
site interactions the (say) 3rd best vector may jump from 
the eigenvector number 3 to the eigenvector number 4, a 
problem which will be documented below. In such a case 
the effective Hamiltonian will be a discontinuous function 



®*abc = ^abc^abc an d from the spectral definition of 

H e?r 

H L ABC = E*ABc\^ABc)(^ABc\+E*ABc\^*ABc)(^iBc\ 
+E*ABC\<B C ){<BCI (36) 

one may calculate the diagonal matrix elements of H ABC 
and reexpress its matrix elements as 

(r A ^c\H e A f i c \rA^c) =e a +e%+e° c + 

V(A*)B + V(A*)C + VBC +V( A *)BC, (37) 

which defines a three-body interaction V(a*)bg, an d re- 
vised excitation hopping integrals 

Wa^c^Mc^Pb^c) - h AB + h AB{C y (38) 

The last term represents the effect of C on the hopping 
between A and B. One also obtains effective hopping 
between non directly interacting blocks (for instance A 
and C through B in a linear ABC configuration). This 
indirect propagation may proceed, for instance for triplet 
states, through the process V'aV'bV'c < — > ^a^b^c < — y 
ip A ip B ^Q. These effective interactions are used in the ex- 
citonic treatment. 

For the periodic lattice 



(39) 
(40) 

I 

of the intersite interactions, which is a rather unpleasant 
feature. 



D. Comment 

Of course the method is anly applicable to the study 
of gapped systems and to locate the phase transition be- 
tween a gapped phase and a gapless phase. This limit is 
due to the fact that one uses different model spaces for 
the ground state and for the excited states. The method 
cannot provide the low energy physics of gapless anti- 
ferrmagnetic lattices. For such phases the method is un- 
able to give a strictly zero gap nor the density of states. 
As will be shown in the following examples the calculated 
gap becomes extremely small and in some cases it may 
be spuriously negative. This limit (which is not present 



I 

= (■■■^u^j---\H eff \--^u^j---) = J2 E "+ E i+ E E v ^ 

J±I J&I) <K&I) 

E E E v j kl+ E E v (i-)jk, 

= (■■■ 1> hWj ■ ■ ■ \H eff I ' ' ' tPUiPj ■■■) = hu + h u(K) ■ 

K 
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in CORE method) should be kept in mind. 



III. TEST APPLICATIONS 

A. Illustration of the superiority of effective 
interactions over bare ones 

This first subsection will illustrate the crucial effect of 
considering effective interactions in excitonic models. We 
shall consider textbook problems concerning the 1-D non- 
dimerized chain using a tight-binding mono-electronic 
Hamiltonian. 



H = t 2J (cip a q + a+a p ). 



(42) 



We shall address successively three problems, namely, 

- the position of the highest orbital for half filling, 

- the energy of the lowest occupied orbital, and 

- the size-dependence of the lowest excitation energy 
at half filling. 

In the three problems we shall demonstrate that a bare 
excitonic model gives an incorrect behavior of the prop- 
erties, introducing spurious dependence on the size of the 
blocks. On the contrary the definition of effective inter- 
actions between the blocks dramatically reduces the error 
coming from the use of finite size blocks. 

One may address first the question of the Fermi level 
for the half-filled band, which is of zero energy. One shall 
consider blocks of size 2n. In each of them (say I) the 
ground state is 



$ = n # 



(43) 



fe=l,- 



where ffj^f = £ k (f> k I with £ k = 2t cos 

For the whole system the ground state is treated from 



*°=rw- 



(44) 



The lowest ionization process in each block concerns the 
highest occupied orbital <fi™, 

4>i = a^°j. (45) 

The lowest ionization energy for the block is 



AJ = -£ r ' 



-2t cos ■ 



2n + 1 2n + 1 



0(2), (46) 



where 0(2) is proportional to n~ 2 . The ionized state of 
the whole system will be described according to Eq. fTTI 
from the vectors ^} = a^^°. The direct matrix element 
of H between ^} and is 



b ^°\H\a^°) = -(M\H\rj) 



(47) 



Since 



(48) 



where \\ is the orbital localized on site r in block /, 



and since 



- -C n n C n (xi n \H\xj/ 

1\2. 



™. 1 . nnr 

CL = , sm ■ 



" V^TT 2n + 1 



(49) 



(50) 



t . 717T 2 t 



(1 + 0(2)). (51) 



This is the bare interaction, propagating the hole between 
adjacent blocks. The bare excitonic method predicts 
therefore the ionization potential of the infinite chain cal- 
culated blocks of 2n sites as 



A*°°(2n) = A* I + 2h IJ 
tir 



2t 



2n + l 2n + l 



f(-7T + 2) 

2n+ 1 



(52) 



This energy is smaller than that of the 2n sites block but 



it remains different from zero and of order 



On the 



contrary, the effective interaction can be calculated from 
the energies of the highest occupied orbitals of the An 
sites system. It will lead to a specification of Eq. E21 

A°°(n) = 2A*(4n) - A*(2n) 

The spurious n _1 dependent term, present in the bare 
excitonic model, disappears from the effective one, as de- 
sired. 

As another application one might have looked at the 
lowest energy level, which should be 2t for the infinite 
lattice. For a 2n sites block the corresponding energy is 



£ 1 (n) = 2£cos 



2n+ 1 



(54) 



The deviation to the asymptotic value is n 2 t/4n 2 . The 
bare interaction 



t 



t 



2n+l 



2n + 1 



(55) 



behaves as n~ 3 and cannot compensate the n~ 2 de- 
pendent error which affects the £ 1 (n) value. On the 
contrary the effective interaction between these two 
blocks is (l/2)(£ 1 (4?i) - £ 2 {An)) ~ 3iT 2 t/32n 2 i.e., has 
the expected n~ 2 dependency and reduces the error to 
ir 2 t/32n 2 . This problem has an historical importance 
since the recognition of the failure of the original RSRG 
version on this problem has led to the development of the 
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2 nd order 



intra-bloct inter-blocli intra-block 
excitations excitations excitations 



4 th order 



higher order 

FIG. 1: Schematic view of the incorporation of complex effects 
in the off-diagonal effective interactions. 

density matrix renormalization group (DMRG)iA&*A£ One 
sees that the consideration of properly extracted effective 
interactions represents a qualitative improvement of the 
RSRG method, as already mentioned elsewhere 

As an even more dramatic problem, one may con- 
sider the excitation gap in the same 1-D problem, at 
half-filling. In a 2n sites block the local excitations are 
(<b n — > <b n+1 ) electron jump around the Fermi level, the 
corresponding excitation energy is A} = — 2t/{2n + 1). 
There is no bare interaction between the state which is 
excited on block I and on the one excited on block J, 
due to the mono-electronic nature of H. 



b n^ \H\a+ 1)fl ^$o) =0. 



(56) 



Hence the excitation energy would be that of an isolated 
block and would behave asn -1 . On the contrary there is 
a n^ 1 dependent effective interaction which cancels the 
spurious n" 1 term since Eq. 133 leads to 

2t?: 4i?r 



The term 



2n + 1 An + 1 



+ 0(2) = 0(2). (57) 



+ 



effl 



f*0> 



2tn 



(58) 



(4n+ 1)' 

reflects the indirect interaction through inter-blocks 
charge transfer states, in particular of the type 
a + {n+1) a^^>o, which make possible the independent de- 

localization of the hole and of the particle (cf. Fig. 

B. The 1-D frustrated spin chain 

The 1-D spin anti-ferromagnetic (AF) chain with J\ 
spin couplings between nearest neighbor sites and J 2 cou- 



plings (also AF) between next-nearest neighbor sites (cf. 
Fig. |2j, is ruled by the Heisenberg Hamiltonian 



H 



2Ji S i S i + i + 2J 2 S i S j + 2- 



(59) 



It presents a second order phase transition for (J 2 / ' Ji) c = 
j c = 0.2411i 1 7i 1 ff. 1 i? There is no gap for J 2 /Ji = j < j c 
while a finite gap exists beyond this critical ratio. Close 
to the critical point the gap increases very slowly, pre- 
senting as essential singularity at j c . It behaves at this 
origin^ as 



(3 exp • 



] - Jc 



(60) 



DMRG calculations have been reported for this system^ 
as well as analytic treatments^ The renormalized exci- 
tonic method has been applied to (n = 4, 6, 8, and 10 
sites) blocks, and extrapolated. For a given value of n, 
the calculated gap A*°° for the lattice, estimated from 
Eq. EH! using the A*(n) and A*(2n) excitation energies, 
is dramatically reduced with respect to A*(2n), due to 
the cancellation of the (n -1 ) components of the A*°° ex- 
citation. Actually in such a simple problem 
A*°° = 2A*(2n)- A*(n). If A*(n) = A + Bn^ + Cnr 2 



'{n)=A + C{ 



1 



1 



4n 2 16n 2 ' 



A 



3C 
16n 2 



(61) 



From the different calculations of A*°°(n) it is possible 
to estimate an extrapolated value of the gap. We have 
used a polynomial fit 

A*°°(n) = ai + a 2 (n + 1)~ 2 + a 3 (n + I)- 3 + a 4 (n+ l)" 4 , 

(62) 

which gives the results reported in Fig. 03 One may 
notice that 

- the extrapolated value of the gap A*°° = a\ for 
j < j c is not strictly zero. The largest error is 
for j = where A*°° = 0.0068 J x . This value is 
within the accuracy of the extrapolation techniques 
of DMRG (cf. Fig. Elof Ref. 19 ), 

- the calculated gap goes through a minimum at 
j = 0.24, close to the critical value, where |A*°°| = 
3.10~ 6 

- it increases for larger values of J 2 . Immediately 
beyond J 2c the gap follows the expected law. We 
found 

a = 0.21022 
(3 = 0.10253 




FIG. 2: The non-dimerized frustrated 1-D chain. 
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FIG. 3: Dependence of the gap of the frustrated 1-D chain 

on the j — J2/J1 ratio. o and □ direct 

gaps from 16 and 20 sites segments, • • • o • • • and 
REM gaps from 8 and 10 sites blocks. The arrows indicate 
the benefit of the REM treatment. The full line gives the 
extrapolated gap from REM. 



for the parameters of Eq. Efll 

- the calculated gap for the Majumdar- Ghosh point 
(2J2 = Ji) is 0.465Ji, which compares well with 
the DMRG^ estimate (0.48 Ji) and the result of 
an analytic development (~ 0.45Ji)2fli 



C. The 1/5-depleted square 2-D lattice 

The 1/5-depleted square 2-D lattice, built from square 
plaquettes and octagons (cf. Fig. 0J, was first considered 
as representing the 2-D lattice of the CaViOg crystal. It 
appeared later on that next-nearest neighbor spin cou- 
plings are important in this material, but the simple pic- 
ture, with J p AF couplings on plaquette bonds and Jd AF 
couplings between adjacent plaquettes already presents 
an interesting physics with three phases. When the pla- 
quettes are weakly coupled, i.e., J P /{J P + Jd) = j > j c , 
the system, in this plaquette phase, is gapped. It is also 
gapped when the dimers connecting the plaquettes are 
weakly coupled, i.e., when j < j' c . This phase is called 
dimer-phase. In between, i.e., for j c < J P /(J P + Jd) < j' c , 
the lattice keeps a Neel order and this phase is gapless. 
Several studies, using perturbative expansions^! or quan- 
tum Monte Carlo (QMC) calculations^ agree on this pic- 
ture and propose j c ~ 0.4±0.01 and j' c ~ 0.51±0.01. We 
have tested our method on this problem. The simplest 
block that one may consider is the octagon (see schema 
(A) of Fig. 0}. It is non-degenerate whatever the J p j Jd 
ratio. Actually starting from these blocks, REM provides 
a correct picture of the physics, since the gap disappears 
between j c = 0.40967 and j' c = 0.50945 (cf. Fig. EJ. This 
result is obtained from 8 sites blocks. 
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FIG. 4: 1/5-depleted 2-D square lattice, definition of various 
blocks. 



In order to ckeck wether this excellent agreement was 
not fortuitous we have introduced next-nearest neighbor 
interactions between octagons, applying the formalism of 
section II C. Two types of trimers (linear and perpendic- 
ular) have to be considered. The results appear in Table 
1, and they deserve the following comments 

- the dependence of the gap on the j ratio is almost 
the same as when working with dimers only. The 
gapless domain in slightly reduced to the interval 
0.39572 < j < 0.49784 

- the 3rd target vector for the perpendicular trimer 
(i.e., th 3rd vector presenting the largest projection 
on the model space) is the 3rd eigenvector (^3) 
of the perpendicular trimer problem for j < 0.41 
and the 4th (^4) one for j > 0.41. This may be 
seen as a signature for a finite (24 sites) cluster of 
the vicinity of the phase transition in the periodic 
lattice. A similar phenomenon is observed for the 
linear trimer between j = 0.40 and j = 0.41. Re- 
garding the Neel-plaquette phase transition, a sim- 
ilar change of the target vectors appears for 0.50 < 
j < 0.51 in the perpendicular trimer superblock. 
This phenomenon of discontinuity of H e f* does 
not appear when working with dimers only. One 
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TABLE I: Calculated spin gap for the 1/5-depleted square 
lattice as a function of j = J v j (J p + Jd) 
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might eventually circumvent this problem by tak- 
ing a weighted energy for the third "root" appear- 
ing in the spectral definition of H e $' . If Pq is 
te projector on the model space a = ||Po*3|| = 
(P * 3 |Po*3), /3 = \\Po*4\\ = (fb*4|Po*4> and 
if H AB c\^3) =E 3 V 3 ), H ABC \t>4) = E 4 V 4 ), one 
might define E* ABC — aEs + f3En, i.e. proceed to a 
diabatization of the 3rd target vector. 



One may alternatively change the shape of the blocks. A 
stared 8 sites block with 4 dimer bonds around a plaque- 
tte pictured in scheme (B) of Fig. 01 is expected to be 
relevant for the dimer phase. The gap calculated from 
these blocks almost coincides with the previously calcu- 
lated one, with a critical value j c = 0.40 (see Fig. EJ, 

For the plaquette phase we have considered blocks with 
one or two plaquettes, all atoms belonging to plaquette 
bonds. The number of sites are 8 and 10 respectively, as 
pictured in schemes (C ) and (C) of Fig. 0] An extrapo- 
lation leads to the curve on the right part of Fig. The 
gap is slightly smaller than from the octagons, but the 
critical values of disappearance of the gap, j' c — 0.51, co- 
incide. We have compared our calculated gaps with the 
ones reported (in Figure 2 of Ref. 22 ) from QMC calcu- 
lations, and the two methods practically coincide (within 
the uncertainties of reading of the above-mentioned Fig- 
ure). 




Jp+Jj 



FIG. 5: Gap in the 1/5-depleted 2-D square lattice. ( — ) from 
octagonal blocks (A), (...) from blocks of type (B), (- -) from 
blocks of type (C') and (C), after extrapolation^ o) QMC 
calculations from Ref. 22. 



D. The plaquette lattice 

The square type lattice built from interacting pha- 
quettes is characterized by intra-plaquette J and inter- 
plaquette j AF couplings (cf. Fig. HjJ. The properties de- 
pend on the A = j/J ratio. For j = 0, the plaquettes are 
independent and the lattice is gapped. It is not gapped 
for the j = J 2-D square lattice, and phase transitions as 
expected to occur for (j/J) c = A c and for (J/j) c = 1/A C 
(this last relation being due to the intrinsic symmetry 
between j and J). Several works have been devoted to 
this problem. Third order series expansions^ and QMC 
calculations^ suggest that A c ~ 0.55. Extrapolations of 
finite size exact diagonalizations^ fail to give a zero spin- 
gap whatever the value of A. A recent work has used the 
CORE method^ together with order parameter suscep- 
tibilities, suggesting a critical behavior between A = 0.5 
and A = 0.6. The problem of the gap is reexamined here 
using REM. 

Two types of blocks have been considerd. 
The first one involves one, two or three plaquettes (n = 
4, 8 and 12 sites), fragments of a ladder (see schema (A) 
of Fig. Hjjl, There are two types of dimers, collinear or 
side by side. If one calls AEab and AE' AB the excitation 
energies for these dimers, AE A the excitation energy of 
the block, our model leads to the following expression 
of the gap AE(m,n2), for a block of n\ sites along the 
longitudinal direction, 71% sites along the transverse one, 

AE(m,n 2 ) = 2AE AB (2n u n 2 ) +2AE AB (n 1 ,2n 2 ) 

-3AE A (m,n 2 ). (63) 

Fig. reports the gap calculated for n x — 2,3,4 and 
n 2 = 2. One sees that the gap vanishes when A tends to 
1. An extrapolation is possible in terms of nj -1 and nj -2 , 



FIG. 6: Definitions of blocks for the study of the plaquette 
lattice. 



for U 2 



AE(ni, n 2 = cte) = Aq 



Ai 

Tlx 



(64) 



AE(n 2 ) = A is the ni extrapolated value of the gap for 
a fixed value of n 2 . Assuming that 



7 

AE(m,n 2 ) = (a+ — + -^- 
ni nf 



_ a-i a-r a 2 b 

AE{n u n 2 ) = AE + — + — + -§ + -4 



n\ n 2 n\ n 2 n\n 2 



one obtains 



AE(m,n 2 ) = AE "- 



2n\ 2n 2 n\n 2 



(66) 



(67) 



Confronting Eq. Eland Eq. IfTTI gives A = AE - a 2 j2n\, 
A\ = —b/n 2 and A 2 — —a 2 /2. Hence the final value 
of the extrapolation gap is AE = A Q - A 2 jn\. Fig. 
reports the so-calculated gap as well as the values 
AE(ni,n 2 ) for n 2 = 2 and n\ = 2,4,6 as a function 
of the j/(J + j) ratio. One sees that a gapless phase 
appears for n\ > 2. After extrapolation the lattice is 
found to be gapless for j/J > 2/3 = 0.666. The value of 
A c is somewhat larger than the commonly accepted value 
but it represents a considerable improvement over the ex- 
trapolations of finite lattices. One may mention that, as 
a by-product of the present calculations, one obtains, for 
j = J, a value of the gap of the two-leg ladder. The ex- 
trapolation leads to AE = 0.47J, close to the best QMC 
estimate (0.50J)i2i 

A second type of rectangular blocks have been consid- 
ered, involving odd numbers of sites in one direction and 
even numbers in the other one (see schema (B) of Fig. 
El- In such a case there are three types of dimers. These 
12 sites (ni = 3, n 2 = 4) blocks are more compact than 
the previous (rii = 6, n 2 = 2) ones and the calculated 
gap, which appears in Fig. is somewhat lower. Ex- 
trapolation is difficult in this case, due to the difference 
in the physical nature of the dimers, but the evaluations 
from different blocks are quite consistent. 



J 

J+j 



FIG. 7: Singlet-triplet spin gap in the plaquette lattice from 
2x2 block (A), 4x2 block (□), 6x2 block (+), 4x3 block 
(o). Dashed line: extrapolation from n x 2 block. 



IV. CONCLUSION 

We have presented a very simple method for the study 
of the gap in gapped periodic lattices. The method rests 
on the consideration of blocks and a truncation of the 
Hilbert space to products of a few eigenstates of the block 
as practiced in the RSRG. In the past we have considered 
(2n + 1) sites blocks, with spatially non-degenerate dou- 
blet ground states, in spin lattices. The blocks can then 
be seen as S z = ±1/2 quasi-spins. Using the theory of ef- 
fective Hamiltonians, and the exact spectrum of dimers 
(or trimers) of blocks, we have proposed to renormal- 
ize the interactions between blocks, and the so-obtained 
variant of CORE (RSRG-EI) happens to keep the con- 
ceptual elegance of Wilson's idea while gaining, at a very 
low cost, numerical relevance 

The present work is closely related but different. It 
considers blocks with even number of sites, presenting a 
non-degenerate ground state. Again the exact treatment 
of the block and of the dimers or trimers of blocks is em- 
ployed to define block effective energies and inter-block 
effective interactions. However different model spaces 
are used for the ground state and for the lowest excited 
states. For the ground state ^o, built from the product 
of ground states, the energy is a simple sum of intra and 
inter-block energies. The excited states are linear combi- 
nations of locally singly excited functions, products of an 
excited state on one block by the ground states functions 
on the other blocks. This space is a small fraction of those 
handled in RSRG techniques. The knowledge of the ex- 
cited states of dimers or trimers of blocks enables one to 
define the effective interactions between an excited block 
and its ground state neighbors, as well as effective excita- 
tion hopping integrals, which delocalize the excitations. 
The effective interactions incorporate complex processes, 
including multiple excitations or/and inter-block excita- 
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tions. These informations are used to build an excitonic 
Hamiltonian for the infinite lattice, and to estimate the 
gap. 

The method does not provide any information on the 
low energy physics of gapless phases but solves some dra- 
matic failures of the use of bare interactions (as they 
manifest in the original RSRG formalism). It can be ap- 
plied to various Hamiltonians (Tight-binding, Hubbard, 
Heisenberg, ab-initio). The method has been presented 
(and tested) in its simplest version on spin lattices, with 
identical blocks, one excited state per block, and extrac- 
tion from dimers and trimers. It is possible to generalize 
it to blocks of different sizes or topologies, and one may 
keep several excited states per block. The renormalized 
excitonic method has been tested so far to the research of 
singlet-triplet gaps but it is applicable as well to singlet 
to singlet excitations. The bottleneck is the size of the 



dimers or trimers of blocks, the lowest states of which 
have to be calculated. In the few benchmark problems 
tested in the present work the results are surprisingly ac- 
curate and the method seems to be able to locate phase 
transitions between gapped and gapless phases in 1-D 
and 2-D lattices at a very low computational cost. An 
other application concerning the Shastry-Sutherlandr& 
shows the relevance the here proposed method for the 
study of phase transition in frustrated 2-D spin lattices. 
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